Fractal Structure of Hastings— Levitov Patterns Restricted in a Sector Geometry 
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A generalized form of the Hastings and Levitov (HL) algorithm for simulation of diffusion- limited 
aggregation (DLA) restricted in a sector geometry is studied. It is found that this generalization 
with uniform measure produces "wedge- like" fractal patterns in the physical space, whose fractal 
dimension and anisotropy exponent depend significantly on the opening angle f5 of the sector. The 
morphological properties and the overall shape of the patterns are analyzed by computing the angular 
two-point density correlation function of the patterns. We also find that the fractal dimension of the 
patterns with sinusoidal distributed measure depend weakly on f5 with almost the same dimension as 
the radial DLA cluster. The anisotropy exponent and the visual appearance of the patterns in this 
case are shown to be compatible with those of the advection-diffusion-limited aggregation (ADLA) 
clusters. 

PACS numbers: 61.43.Hv, 82.40. Ck, 68.43.Jk, 64.60. al 
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I. INTRODUCTION 

Diffusion-limited aggregation (DLA) was originally in- 
troduced by Witten and Sander [1] in 1981 to model the 
aggregates of metal particles formed by adhesive contact 
in low concentration limit. This model has been then 
shown to describe many pattern forming processes in- 
cluding dielectric breakdown [2], electrochemical deposi- 
tion [3, 4], viscous fingering and Laplacian growth [5]. 
One of the standard approaches to simulate a Laplacian 
field is by random walkers, which are launched from the 
periphery of the system and diffuse toward the growing 
cluster and freeze on it. This procedure is equivalent to 
solving Laplace equation outside the aggregated cluster 
with an appropriate boundary conditions. The walker 
sticks to a point on the surface of the aggregate with a 
probability proportional to the harmonic measure there 
[6]. 

Another powerful method for studying such growth pro- 
cesses in two-dimensions, is the iterated stochastic con- 
formal mapping [7-9], which is known as Hastings and 
Levitov (HL) method. 

Since DLA enhances the instability of growth, the re- 
sulting cluster is highly ramified and branched. The pro- 
cedure of the proliferation of the fingers in the fractal 
structure of DLA is one of its complexity sources which 
is one of the purposes of the present paper to deal with. 

The HL method is based on the fact that there exists 
a conformal map that maps the exterior of the unit cir- 
cle in the mathematical i(;-plane to exterior of the cluster 
of n particles in the physical z-plane. The complicated 
harmonic measure in the physical space then changes to 
a simple measure: a uniform distribution in the mathe- 
matical space. The HL algorithm can be also used for 
more generalized growth models for which the probabil- 



*Electronic address: a_saberi@ipm.ir 



ity measure is not uniform in the mathematical plane. 
In a part of this paper, we also study the statistical and 
morphological properties of the fractal patterns produced 
by a novel generalization of the HL algorithm with both 
uniform and non-uniform probability measures in the w- 
plane. We restrict the algorithm to attribute the mea- 
sure in the i^-plane to a sector of opening angle /3. We 
find that this restriction which breaks the radial symme- 
try in the growth process, yields fractal patterns whose 
dimensions depend on whether uniform or non-uniform 
distribution is used for the measure. For the uniform dis- 
tribution, the patterns looks like those of DLA in a wedge 
but with fractal dimension depending on the opening an- 
gle p. However for the non- uniform distribution of the 
measure, we find that the statistical properties of the pat- 
terns is compatible with advection-diffusion-limited ag- 
gregation (ADLA) [10] and their fractal dimensions take 
values close to that of radial DLA i.e., df = 1.71 with a 
weak dependence on the opening angle /3. 

Bazant et al. [10] proposed a method based on time- 
dependent conformal maps to model a class of non- 
Laplacian growth processes such as ADLA in a back- 
ground potential flow. They showed that in spite of dra- 
matic increases in anisotropy, the fractal dimension of 
the patterns is not affected by advection and takes the 
same value as for radial DLA clusters. 

The fractal dimension of the patterns is already shown 
to be affected by the geometrical factors. Stepanov and 
Levitov [11] obtained dimensions as low as df = 1.5 for 
simulations of anisotropic growth (a different anisotropy 
that we consider in this paper) using the noise-controlled 
HL algorithm. A similar result reported in [12] for DLA 
with anisotropic perturbations. Davidovitch et al, [13] 
introduced growth models that are produced by deter- 
ministic itineraries of iterated conformal maps with dif- 
ferent dimensions. In channel geometry, the fractal di- 
mensions df = 1.67 [14] and df = 1.71 [15] are reported 
for periodic and reflecting boundary conditions, respec- 
tively. 



The comparison of the ensemble average of DLA cluster 
density with the noise-free Saff man- Taylor viscous fin- 
gers [5] have been the subject of various studies in both 
channel and wedge geometries [16-19]. 

There has been shown in [20] that there exists a critical 
angle r] ~ 60° — 70° in viscous fingers and DLA growing 
in a wedge, indicative of a typical angular spread of a 
major finger. In this paper, for the patterns produced by 
restricted HL algorithm with uniform measure, we find 
that their visual appearance like a DLA cluster in a wedge 
and we compute the average angle between two major 
fingers. To estimate the angle, we measure the angular 
density-density correlation function. We also obtain the 
relation between the minimum opening angle of a wedge 
which contains the cluster in the z-plane, and the opening 
angle /3 considered in the w-pldine. 



II. RESTRICTED HL ALGORITHM IN A 
SECTOR GEOMETRY 

Hastings and Levitov (HL) employed the conformal- 
mapping tool to describe the complicated boundary of a 
growing DLA cluster. Their approach [7] was based on 
iteratively applying the function (/)a,6>('^) which maps a 
unit circle to a circle with a bump of linear size vA at 
the point w = e^^ ^ 
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The parameter < a < 1 determines the shape of the 
bump, for higher a the bump becomes elongated in the 
normal direction to the boundary dC (e.g., it is a line 
segment for a = 1). In this paper we set a = ^ for which 
the bump has a semi-circle shape. 

A cluster Cn consisting of n bumps can be obtained by 
using the following map on a unit circle 

which corresponds to the following recursive relation for 
a cluster Cn+i, 



^n+l(^) =^ 
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In order to have bumps of fixed-size on the boundary 
of the cluster, since the linear dimension at point w is 
proportional to |<I>^(il')|~^, one obtains 
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To produce an isotropic radial DLA, the harmonic mea- 
sure of the n-th growth probability p(z, n) has to be con- 
formally mapped onto a constant measure p{6n) = ^ 



on a unit circle i.e., < On < ^^tt (more details of the 
simulation algorithm which have been also considered in 
this paper are given in our previous work [21]). 

In this paper we use a generalized form of the HL al- 
gorithm in which ^^s are restricted to be selected from 
a fixed interval < On ^ P on di unit circle, with both 
uniform and non-uniform measures. This measure 
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otherwise, 
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induces a complicated harmonic measure on the bound- 
ary of the patterns in the 2;-plane. 



HI. SCALING PROPERTIES OF THE 
PRODUCED FRACTAL PATTERNS 

The function ^n{uj) in Eq. (3), has the following Lau- 
rent series expansion 
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where F| ^ is called the conformal radius and Fq ^ the 
center of charge of the cluster Cn- It is possible to analyt- 
ically show that [9] these coefficients contain descriptive 
information about the morphology of the clusters. In 
particular, for the coefficient of the linear term in Eq. 
(7), it can be shown that [9] 
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which scales with the cluster size n as 
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where df denotes the fractal dimension of the cluster. 

Scaling behavior of the next Laurent coefficient F^ 
can also explain about the isotropicity of the clusters. For 
an isotropic DLA cluster, this coefficient scales according 
to the following relation 
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with 2/(io = 0.7 [9]. This also holds for the anisotropic 
growth phenomena but with a different exponent do, e.g., 
do - 1.71 is obtained for ADLA [10]. 
Since our produced patterns are anisotropic, in order to 
have a quantitative measure and compare them with that 
of the other studied models, we check the scaling rela- 
tion (10), and compute the exponent do as a function of 
the opening angle P (for convenience, we address do as 
anisotropy exponent^ in this paper). We compute Fq^ by 
using the recursion equation given in [9] for the Laurent 
coefficients. 



IV. RESTRICTED HL PATTERNS WITH 
UNIFORM MEASURE 



In this section, we investigate the effect of the resti 
tion made on the angular distribution of the measure 
the mathematical plane, on the statistical properties 
the patterns produced by using the algorithm descrit 
in section II, with p{9n) = ^ and < 6n < /3. 
We grew 400 clusters of size n = 10^, and 30 clust 
of size n = 10^ for different opening angles in the rai 
4.5° < /3 < 360°. Having looked at the sample pattei 
shown in Fig. 1, the appearance of a wedge- like shg 
in the all patterns is evident. Each of the two inter 
limits for the distribution of ^^s, i.e.. On = and /3, can 
represent itself as a boundary of the wedge (which due to 
the conformal maps, are not necessarily straight lines) in 
the physical plane. Despite that the harmonic measure 
in iL'-plane is uniform, the growth probability looks much 
larger around the boundaries in the 2:-plane. For small 
angles of /3 < 36° there always appears only one finger 
with sidebranches on a boundary. For larger /3 there also 
exists at least one finger which has selected a boundary 
to grow on. In the range angle that two main branches 
coexist, because of the attraction of the boundaries, the 
average angle between branches increases by increasing 

The snapshots of an example of growing cluster with 
f3 = 216° in Fig. 2 graphically show the procedure of 
the branch formation in the growing fractal pattern. 



(a) 






FIG. 1: Some typical clusters of size n — 10^, generated by 
using the restricted HL algorithm with uniform measure and 

(a) /3 = 36°, (b) P = 324°, (c) /3 = 144° and (d) /3 = 216°. 



> 





n = \(f 



n = 3xl0^ 



n = 5xlO^ 



FIG. 2: Snapshots of a growing cluster of (5 ■ 
number of iterative conformal maps n. 
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A. Scaling exponents 

To have a quantitative understanding of the behavior 
of our clusters, let us now measure some of their sta- 
tistical quantities. The first quantity we would like to 
measure is the fractal dimension which can be obtained 
using the scaling relation in Eq. (9). Two sets of re- 
sults, one averaged over 400 clusters of size 10^ and the 
other over 30 clusters of size 10^ at each opening angle 
/3, are reported in Fig. 3. The fractal dimension shows a 
significant dependence on the opening angle /3. For nar- 
row angles the dimension takes values very close to unity 
which is dominated by the boundary effects. By increas- 
ing /3, it reaches a relative maximum at around /3 c^ 18° 
which is almost the half-maximum angle interval in which 
there exists only one main branch. There also exists a 
local minimum of the dimension around j3 :^ 144° which 
seems to be a typical angle up to which only two main 
branches coexist. For greater angles, the dimension in- 
creases until for /3 = 360° where an isotropic radial DLA 
with df = 1.71 is expected. For DLA clusters in a wedge 
produced by the original definition [1], the fractal dimen- 
sion is not affected by the wedge geometry of opening 
angle a [20], and it depends weakly, if at all, on a. 

The other quantity which we measure for our clusters 
is the exponent describing the scaling behavior of the 
position of the center of charge |Fo | with the cluster size 
n, defined in Eq. (10). We find such a scaling behavior 
only in the interval 9° < /3 < 180°. There has been also 
observed that the corresponding exponent d^ depends on 
P in this range. The results for two mentioned sets of 
averages are shown in Fig 4. 



Overall shape and the morphology of the 
clusters 



In this subsection we discuss about the overall shape 
and the morphology of the clusters. As mentioned before, 
for sectors of angle j3 < 36° there exists only one main 
branch, and for greater angles up to /3 ::^ 144° there are 
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n = 10,000 
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density-density correlation function is then read off from 

c{<j>) = Me + <j>)p{e))-{pf]\e, (11) 

where p{9) is the density of particles in the cluster in a 
50 = 1° sector around 6 (see Fig. 5), and (••) denotes 
for averaging over the sectors in 2:-plane. The computed 
correlation functions for different opening angle /3 aver- 
aged over 30 realizations are shown in Fig. 6. As can 
be seen in the figure, for clusters consisting of one main 
branch there is an anticorrelation between the origin and 
other angles. Appearance of the second peak in the func- 
tion with positive correlation is indicative of the second 
coexistent main branch in the cluster. 

Location of the second peak shows approximately the 
angle between two main coexisting branches rj in the 
physical plane which is reported in the second column 
of Table I. 



FIG. 3: The fractal dimension df of the patterns produced 
by using the restricted HL algorithm with uniform measure 
as a function of the opening angle /3 in K;-plane. The averages 
were taken over 400 clusters of size n = 10^ (circles) and 30 
clusters of size n = 10^ (squares). 
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FIG. 4: The exponent do as a function of /3 in the scaling 
region, according to Eq. (10). 



two main coexistent branches, and for angles /3 > 144° 
more than two main branches appear. To quantify this 
observation and obtain the average angle between two 
main branches of the patterns, we follow the same rea- 
soning as in [20], for our clusters. Kessler et al.^ [20] 
studied the building block of DLA clusters in a wedge 
by computing the angular two-point density correlation 
function in the constitutive sectors of a cluster. 

Consider two sectors separated by angle (j). The 
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FIG. 5: Schematic description of the procedure expressed in 
the text to compute the angular correlation function accord- 
ing to Eq. (11). 

For a given opening angle /3 in i^-plane, the produced 
patterns have a wedge-like shape in the 2:-plane. To have 
an estimate relation between two angles /3 and the av- 
erage opening angle of the wedge-like patterns in the z- 
plane a, we report some of these corresponding angles in 
the third column of Table I. The averages are taken over 
30 clusters of size n = 10^ for each reported a. 
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a 


36° 


there is only one main branch 


15° 


90° 


31° - 34° 


44° 


144° 


81° -88° 


97° 


198° 


more than two main branches exist 


164° 



TABLE I: The average angle between two main branches 77, 
and the average opening angle of the wedge-like patterns in 
the 2;-plane a for different opening angle /3 in the K;-plane. 
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FIG. 6: Angular correlation function c(0)/c(O) for clusters 
produced by using the restricted HL algorithm with uniform 
measure for different opening angle /3 as a function of 0//3. 
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FIG. 7: Some typical clusters of size n = 10^, generated by 
using the restricted HL algorithm with sinusoidal measure i.e., 

sm{7v0n/l3) distribution and (a) /3 = 90°, (b) /3 = 180°, (c) 
/3 = 270° and (d) /3 = 360°. 



V. RESTRICTED HL PATTERNS WITH 
SINUSOIDAL MEASURE 

As discussed in section IV, the boundary effects are 
evident in the visual appearance and statistical proper- 
ties of the patterns. So this motivated us to force the 
measure to be distributed away from the boundaries i.e., 
On = and /3 in the iL'-plane. By inspiration from [10], we 
produced patterns according to the algorithm described 
in section II, with p{On) ~ s'm{7r6n//^) for < 6n < f^- 
We generated 400 clusters of size n = 10^, and some of 
size n = 10^ for different opening angles in the range 
36° < P < 360°. Some of the clusters are shown in 
Fig. 7. There exists a significant difference between the 
overall shape of the patterns in figures 1 and 7. Due to 
the no n- uniform measure, the visual appearance of the 
patterns in Fig. 7 is not affected by the boundary. 

This difference can be also observed in the fractal di- 
mension of the patterns. As plotted in Fig. 8, except 
for small opening angles, the fractal dimension depends 
weakly on P and takes values close to df c^ 1.7. 

Due to the sinusoidal distribution of the measure, even 
for P = 360° the clusters are anisotropic. The growth 
process is dominated by advection toward a certain di- 
rection with a spatial extent which depends on the open- 
ing angle p. The morphology of these patterns are very 
similar to that of ADLA clusters. 

ADLA is one of the simplest examples of transport- 
limited aggregation (TLA) in which the released random 
walkers are being drifted in the direction of a background 
potential flow. The difference between simulations of 
TLA and DLA by HL algorithm is in the sequences of 
the angles On- In TLA the angles are chosen from a 
time-dependent (non-harmonic) measure p(^,tn). 

We have also examined the scaling relation Eq. (10) 



for our clusters with sinusoidal measure and obtained the 
dependence of the exponent do on the opening angle p. 
As shown in Fig. 8, do is weakly dependent on P with 
values of a little less than df which is roughly in agree- 
ment with the same scaling behavior of ADLA which is 
reported in [10] for a special case p = 360°. 




360 



FIG. 8: The fractal dimension df and the scaling exponent 
do as a function of /3, computed for the patterns generated by 
using the restricted HL algorithm with sinusoidal measure. 



VI. CONCLUSIONS 

The fractal structure and statistical properties of the 
patterns generated by a generalized HL method restricted 
in a sector geometry were studied. It is found that the re- 
striction with uniform measure leads to the production of 
fractal patterns whose fractal dimension and anisotropy 
exponent depend significantly to the opening angle /3 
of the sector in the mathematical i(;-plane. The over- 
all shape of these patterns is governed by a wedge-like 
DLA appearance which has been shown to be the char- 
acteristic feature of the uniform measure. The boundary 
effects lead to a nontrivial dependence of the average an- 



gle between two major coexistence branches on the angle 
/3. The proliferation of the main branches was studied by 
computing the angular density-density correlation func- 
tion giving a quantitative understanding of the morphol- 
ogy of the patterns as a function of p. 
It is also found that the sinusoidal distribution of the 
measure on the sectors gives fractal patterns with al- 
most the same fractal dimension for all values /3. The 
anisotropy exponent and the visual appearance of these 
patterns are shown to behave very similar to those of 
ADLA clusters. 
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